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Abstract 

Lattice models are crucial for studying thermodynamic properties in many physical, biological 
and chemical systems. We investigate Lattice Restricted Primitive Model (LRPM) of electrolytes 
with different discretization parameters in order to understand thermodynamics and the nature 
of phase transitions in the systems with charged particles. A discretization parameter is defined 
as a number of lattice sites that can be occupied by each particle, and it allows to study the 
transition from the discrete picture to the continuum-space description. Explicit analytic and 
numerical calculations are performed using lattice Debye-Hiickel approach, which takes into account 
the formation of dipoles, the dipole-ion interactions and correct lattice Coulomb potentials. The 
gas-liquid phase separation is found at low densities of charged particles for different types of 
lattices. The increase in the discretization parameter lowers the critical temperature and the 
critical density, in agreement with Monte Carlo computer simulations results. In the limit of 
infinitely large discretization our results approach the predictions from the continuum model of 
electrolytes. However, for the very fine discretization, where each particle can only occupy one 
lattice site, the gas-liquid phase transitions are suppressed by order-disorder phase transformations. 







I. INTRODUCTION 



Electrostatic interactions are important in various physical, chemical and biological 
processes.— However, a full thermodynamic description of systems with charged particles 
is still far from complete. In last decade this subject has attracted a lot of attention due to 
controversial theoretical and experimental issues on the nature of phase transitions in ionic 
fluids.^ 

Lattice models have been used extensively for investigations of different phenomena in 
chemistry, physics and biology. For example, the Ising model, which is a lattice gas model, 
is fundamental for understanding critical phenomena in non-charged systems. This obser- 
vation has strongly stimulated many studies of ionic fluids utilizing the discretized lattice 
j^o^glg_M.i2.2ilUU2jL^.i4^ There are several advantages of using lattice approach for sys- 
tems with the particles interacting via long-ranged Coulomb potentials. The production 
of ion pairs and dipole-ion interactions can be described better than in continuum-space 
models.'^'i^ Lattice models are also computationally much more efficient with respect to 
the time and length scales.— i^^i^iii Monte Carlo computer simulations of charged systems 
on discretized lattices are faster by a factor of 5-100 than the corresponding continuum 
models.i^^ 

If we define a diameter of a charged particle as a and the lattice cell size as a, then a 
discretization parameter C, = (j/a specifies how close the lattice system approaches the con- 
tinuum behavior. The case of C = 1 corresponds to the standard lattice model of electrolytes 
where each particle occupies no more than one lattice site. When the discretization parame- 
ter becomes very large the system of charged particles does not feel anymore the underlying 
lattice and its thermodynamic properties become indistinguishable from the properties of 
continuum electrolytes. 

Theoretical investigations of thermodynamics and phase behavior of lattice models of elec- 
trolytes have followed several different approaches.-^ '^'^i^'-'i^'^i^^i^^'^i'^^ The hierarchical reference 
theory^ utilized the renormalization group methods to calculate thermodynamic properties 
of ionic fluids on lattices. Ciach and coworker a^^i^^i^^ used the fleld-theoretical methods to an- 
alyze different properties of lattice electrolytes. However these methods are mathematically 
very complicated and the predictions, as compared with Monte Carlo simulations,^ '^'^ are 
mostly of qualitative nature. A different theoretical method^iiSiiSiii is based on the physically 



1 



more transparent Debye-Hiickel (DH) approach. In this method the free energies of lattice 
electrolytes are obtained by solving the lattice versions of usual Debye-Hiickel equations for 
general dimensions. This method can provide a reasonable thermodynamic description of 
different lattice models of charged particles, such as charge-asymmetric^^ and anisotropic 
lattices,— as compared with computer simulations and other theoretical approaches . ^^i-^^ 
A full thermodynamic investigation of continuum {( = oo) restricted primitive model 
(RPM), which is a system of equal-size hard-sphere ions carrying positive and negative 
charges, using Debye-Hiickel approach has shown that in the system there are gas-liquid 
phase transitions driven by electrostatic interactions.— However, a systematic study of the 
lattice version of RPM^ (with C = 1) indicates that these phase transitions are suppressed 
on simple-cubic (sc) and body-centered cubic (bcc) lattices because of thermodynamically 
more favorable charge-ordering phase transformations. These theoretical observations are in 
agreement with the predictions from Monte Carlo computers simulations.i^*i^ The different 
thermodynamic properties of lattice and continuum models of electrolytes raised the question 
of crossover phase behavior of ionic systems with intermediate lattice discretization. Ciach 
and Stell ,^^i^^i^^ using mean-field arguments and renormalization group methods, discussed 
some features of phase transitions for these lattice models. However, full thermodynamic 
and phase properties of these systems are not investigated yet. The goal of this article is to 
fill this gap. 

In this paper we study the lattice models of charged particles with different discretization 
parameters by applying Debye-Hiickel approach. In addition, we modify the original lattice 
model of electrolytes^ {( = 1) that had several unphysical features. The paper is organized in 
the following way. In Sec. II we provide a general description of lattice models of electrolytes 
using the Debye-Hiickel method. Specific calculations and discussions of thermodynamic 
properties of charged systems for lattices with different discretization parameters are given 
in Sec. HI. Our results and conclusions are summarized in Sec. IV. 

II. GENERAL LATTICE DH THEORY 

Consider the Lattice Restricted Primitive Model (LRPM) of electrolytes in three dimen- 
sions. Note that our method can be easily generalized in d dimensionSj^ however, for clarity 
we prefer to explain the details using three-dimensional lattices. The model consists of 
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equal number of oppositely charged particles, with charges go or —Qo, correspondingly. Due 
to the overall neutrality, the densities of positive and negative ions are equal to each other, 
p+ = p_ = ipi. The particles have a spherical shape with the charge located in the center 
of the sphere. The diameter of the particle is given by a(, where a is a nearest neighbor 
lattice distance and ( is the discretization parameter: see Fig. 1. Thus, in the special case 
C = 1 each ion occupies a single site, while ( = 00 corresponds to a continuum model. All 
particles interact only through Coulombic potential and hard-core exclusion. Then the free 
energy of the system can be written as a sum of two terms, F = F^'^ + F*^'. It allows to 
calculate the thermodynamic properties of the system, such as chemical potentials and the 
pressure, as given bjAiS 

Since every positively charged particle has a corresponding negatively charged ion, we 
can introduce a new hypothetical particles with chemical potential JIi = /!_ + 'p^- Then the 
two-component system is easily mapped into the one-component system with overall density 
of hypothetical particles given by pi. The conditions for phase transitions are specified by 

/^l (T, Pliq) = /il (T, p-aap) , P{T, puq) = P(T, p.,ap) ■ (2) 

Let us take an arbitrary particle and assume that its charge qi is in the center of coor- 
dinates. Then all lattice sites in the system can be divided into three groups with respect 
to the position of this fixed particle as shown in Fig. 1. The first one consists of just one 
lattice cite, r = (0,0,0), where the charge of the particle is located. The second group is 
made of other lattice sites in the exclusion zone where no centers of other particles can be 
found. All lattice sites outside the exclusion zone are in the third group. To determine 
the average electrostatic potential at each lattice site we can derive the linearized lattice 
Poisson-Boltzmann equations following the standard Debye-Hiickel procedures."^"^ For the 
central lattice site we obtain 

4:71 

Ay.(r=(0,0,0)) = --^g„ (3) 

with vq being a unit cell volume. In this equation Aip is the lattice Laplacian that defined 
as follows,^ 

= - A ^['^(^ + ^) - ^w]' (4) 

CnCi 
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where a is a nearest-neighbor vector and the summation runs over all cq nearest neighbors. 

The number of the lattice sites in the exclusion zone (the second group) is determined by 
the discretization parameter ( and the type of the lattice. For example, for sc lattices with 
C = 1, 2, 3, 4 or 5 these numbers are 0, 26, 92, 250 and 484, correspondingly.^ Since there 
are no charges in the exclusion zone outside of the central lattice site, the linearized lattice 
Poisson-Boltzmann equations have a simple form^^ 

A(^(re.) = 0, (5) 

where r^^ specifies the lattice sites in the exclusion zone but not the central point. 

Similarly, for the lattice sites in the third group, i.e., outside of the exclusion zone, the 
linearized lattice Poisson-Boltzmann equations are given b}'^^ 



A(p{Tout) = /«V(r), (6) 

where Tout describes the lattice points outside the exclusion zone, and = 47r/?pigo/-D is 
the inverse squared Debye screening length with j3 = l/fc^T. 

Linearized lattice Poisson-Boltzmann Eqs. Q, © and © can be written in a unified 

way as 

Aifir) = /tV(r) - AoS{r) - ^ A,S{r - r,,ex), (7) 

where the coefficients Aq and Ai are determined exphcitly from Eqs. Q, © and (0), and 
the summation is over all lattice sites from the second group. Note that 6{r) is a lattice 
delta-function defined as 

{1, r = rn; 
(8) 
0, r ^ ro. 

Eq. (jZj) can be solved via Fourier transformation and it yields the average electrostatic po- 
tential ip{r) at every point of the lattice. It allows then to calculate the average electrostatic 
potential at the central site due to all ions except the one fixed at the origin, 

ipi = ipi{r, x^) - ipi{r, = 0), (9) 

where we introduced a dimensionless parameter x = na. It can be shown that the general 
expression for ipi is given by 

^^ = ^G{x'), (10) 
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where the specific functional form of G'(x^) is determined again from Eqs. Q, ® and 

Since the ideal part of free energy of the lattice model of charged particles is known,— we 
only have to calculate the electrostatic contribution. It can be done by utilizing the Debye 
charging procedure,iS*i^ 



All thermodynamic properties and phase behavior of the system can be obtained from the 
free energy through standard calculations.— 

III. LATTICE ELECTROLYTES WITH DIFFERENT DISCRETIZATION PA- 
RAMETERS 

A. Pure DH Theory for sc, bcc and fee Lattiees with C = 1 

Let us discuss in more detail the three-dimensional lattices with the smallest value of 
the discretization parameter, C = 1- This corresponds to the situation when each particle 
occupies one lattice site. Because in this case there are no lattice sites in the exclusion zone, 
the linearized lattice Poisson-Boltzmann equation ^ can be written in the following form. 




(11) 











A^{r) = K^^{r)-AoS{r). 



(12) 



The solution of this equation is found via Fourier transformation. 




(13) 
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where = (27r) ^ J^^ d^k, and we defined a new lattice function J(k) 
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(14) 



The coefficient Aq can be found with the help of Eq. (jHj) , 



A^(r= (0,0,0)) 



/€V(r = (0,0,0))-Ao = A 




6 J (x2 + 6)/6 - J(k) 



(15) 
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Defining the integrated lattice Green's function as 

1 



P z) 



zJ{k) 



we obtain from Eq. p3|) 



^0 



47rgj _ 
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(16) 



(17) 



Tlius the average electrostatic potential at the origin is given by 



4:7iqi a? 



P 



P(l) 



(18) 



from which, after applying the Debye charging procedure as outlined above, the electrostatic 

— el 

free energy density / can be easily calculated. The ideal lattice gas contribution to the 
free energy density is known,^^ 



/ 



id 



^^0 



.In(l-pt), 



(19) 



where p* = pi^o is a reduced density of the free ions. 

Now we can perform explicit calculations for the ionic systems on sc, bcc and fee three- 
dimensional lattices. The lattice functions J(k) are given by^ 



J(k) 



|(cos fci + cos ^2 + cos fca) (sc), 
cos ki cos k2 cos k^ (bcc) , 

\ (cos ki cos k2 + cos k2 cos k-^ + cos ki cos ^3) (/cc) , 



(20) 



with — TT < ki,k2,k3 < 7T. The corresponding integrated lattice Green's functions can be 
evaluated numerically exactly by using the fact that they can be expressed in terms of elliptic 
integrals.iS^ 

The thermodynamic properties and phase behavior of lattice electrolytes on sc, bcc and 
fee lattices are explicitly calculated from Eqs. (H}, (fTSj) and (fT^ . The predicted phase 
coexistence curves are plotted in Fig. 2, while the critical parameters are given in Table 
1. To analyze the effect of different discretizations we also present in Fig. 2 and Table 1 
the phase diagram and critical parameters for the pure DH model of continuum {( = 00) 
hard-core ionic fluids.— It can be seen that the critical temperatures for lattice models 
are approximately 30% higher than the corresponding continuum value. At the same time 
the thermodynamic properties and critical parameters of the lattice pure DH models are 



approaching the continuum model as the number of nearest neighbors is growing from Cq = 6 
(sc lattice ) to cq = 12 (fee lattice), as expected. It also should be noted that in the lattice 
models the predicted phase coexistence curves have a physically reasonable behavior in the 
limit of high density and T 0, in contrast to the original lattice DH theories.^ This is due 
to the fact that the presented theoretical method correctly describes the average potential 
at all densities of charged particles. 



B. Pure DH Theory for sc Lattice with C = \/2, \/3 and 2 

Now we can examine the lattice models of electrolytes with discretization parameters 
C > 1. Thermodynamic properties will be expressed in terms of the reduced density and the 
reduced temperature, namely, 

pl = {aCrpu and T* = (21) 

Consider first the ionic lattice model with ( = \/2. In this case the exclusion zone around 
any charged particle consists of the lattice central site, where the charge is located, and 
6 nearest-neighbor lattice sites where no other charges can be found due to geometrical 
constraints — see Fig. 1. Then the linearized lattice Poisson-Boltzmann equation ((Zj) yields 

Aip{r) = /€V(r) - Ao5{r) - A^J^ ^(^ " ^-n), (22) 

nn 

where a„„ described the position of the lattice sites in the exclusion zone outside of the central 
site. The solution of this equation can be found again using the Fourier transformation, and 
it is given by 

^(r) = Aoj J (^2 + 6)/6_j(k) + y (x2 + 6)/6-J(k)- ^^^^ 



The boundary conditions to determine the coefficients Aq and Ai are found from Eqs. Q 
and ©: 

ft:V(r= (0,0,0)) -Ao, = «V(r = (1, 0, 0)) - Ai. (24) 



Dvo 

From Eqs. (j^^ and we derive 



, [(6 + x^)(-l + |x^(6 + x^)(-l + P(^)] 

° Dvo -(l + a;2)(6 + x2) + a;2(7 + x2)P(^) ' ^ 



^ ^ _4vrg, [|x^(6 + x^)(-l + P(^)] 

' Di;o -(l+x2)(6 + x2)+x2(7 + x2)P(^)' ^ ' 

Thus the average electrostatic potential at the origin due to all other ions, defined in Eq. 

dni), again can be expressed in terms of the integrated lattice Green's function P{z), 



6 \ „ f ^ f 6 



K — P -^r— - ^(1) + 6A' P — — - p(r 

°^x2 + 6 1x2 + 6/ H \x^ + 6' ^' 



(27) 



where A'j = "^^Aj for j = and 1. After Debye charging procedure, as described above 
in Eq. (fTTj). this equation allows us to calculate the electrostatic contribution to thermody- 
namic properties. 

The exact expressions for the ideal free energy contributions for the lattice models with 
> 1 are not known, but for low-densities the expression |T9|l is still approximately valid. 
Then the thermodynamic properties and phase behavior of the system can be calculated in 
a similar way as was done for C = 1 lattices. The critical parameters for the sc lattice model 
of electrolytes with C = are T* = 0.069 and p* = 0.0054. 

For the system of charged particles on the lattice with the discretization parameter ( = 
\/3 the exclusion zone around any arbitrary ion consists of 19 lattice sites: 1 central site, 
6 nearest-neighbor sites, and 12 next-nearest-neighbor sites. For this model the general 
linearized lattice Poisson-Boltzmann equation ((Zj) has the following form 

Aifir) = /€V(r) - Ao6{r) - ^ 6{r - a„„) - A2Y, - a„nn), (28) 

nn nnn 

with Sinnn defining the positions of next-nearest-neighbor lattice sites. As before, we can 
determine the coefficients Aq, Ai and A2 with the help of Eqs. (jS)) and ©, 

. 2,„/ n n\\ A n _ „2,„/„ _ /i n n\N A n _ ,.2, 



«:V(r = (0, 0, 0)) - Ao, = K'ip{r = (1, 0, 0)) -A,, = K'ip{r = (1, 1, 0)) - A2. 
Dvq 

(29) 

Unfortunately, for the lattice models with ( > \pl the electrostatic part of the reduced 
free energy cannot be expressed only in terms of the integrated lattice Green's function. 
However, for explicit numeric evaluations of the thermodynamic properties the following 
equality can be usedi^i^i 

00 

^ _ j^^^ = 3 j dte-^^'li{t)Im{t)In{t), (30) 

k 
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where Ik{t) is a modified Bessel function of the first kind. Then the thermodynamic calcula- 
tions indicate that the critical parameters for the system of charged particles on the lattice 
with C = V3 are given by T* = 0.063 and p* = 0.0047. 

Similar procedure can be performed for the lattice model of electrolytes with the dis- 
cretization parameter ( = 2. The schematic picture for this model is shown in Fig. 1. 
In this case there are 27 lattice sites in the exclusion zone. The general linearized lattice 
Poisson-Boltzmann equation ((7j) can be written as 

Av9(r) = KV(r) -^o5(r) - ^(5(r-a„,„) -A2^(5(r-a„„„) - A3^(5(r-a„„„,). (31) 

nn nnn nnn 

The solution of this equation leads to determination of average electrostatic potential that 
allows to obtain the electrostatic contribution to free energy density. The resulting ther- 
modynamic calculations yield the following values of critical parameters for this model: 
T* = 0.066 and p* = 0.0060. 

Critical parameters of lattice models of ionic systems with different discretizations, ob- 
tained in pure DH theory, are shown in Fig. 3. It can be seen that the increase in the 
discretization parameter ( generally lowers the values of critical temperature and density, 
although the dependence is non-monotonic. This noisy behavior of critical parameters, es- 
pecially for non- integer (, has also been observed in Monte Carlo computer simulations of 
lattice models with large ( > 5.— 

C. Bjerrum Ion Pairing for Lattices with C = 1 

It is known that the pure Debye-Hiickel theory, that takes into account only free ions, 
is not successful for the description of thermodynamic properties of electrolytes at low 
temperatures.— At these conditions the positive and negative particles have a tendency 
to stick together in ion pairs. The framework of the lattice description of ionic systems is 
very convenient for complete analysis of the process of ion pairing, and it avoids the problem 
encountered in the continuum models of electrolytes.— 

Here we consider the pairing as a reversible chemical reaction of the association of pos- 
itive and negative ions and the production of neutral pairs. This process is controlled by 
equilibrium constant K(T) that can be determined from the densities of different species. 
Ion pairs are specified by a density p2 and a chemical potential p2- In the simplest Debye- 
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Hiickel-Bjerrum (DHBj) approximation we neglect the Coulombic interactions between ion 
pairs and free charged particles.i^*i^ 

The condition of chemical equilibrium between the ion pairs and free ions means that 
H2 = fJ'^ + fi- = yUi- Let us introduce thermodynamic activities of the particles in the 
system: 



2z. 



29i ^ 
A? 



Ml 



9o ^ 



= Af 



(32) 

where Aj denotes the de Broglie wavelengths, and A+ = A_ = Ai = A2; the parameters 
6i define the appropriate internal configurational partition functions. "2*^ The chemical equi- 
librium condition can be expressed in the form of the mass-action law, Z2 = \Kz\, which 
yields for equilibrium constant K{T), 



(33) 



With the help of Widom's potential-distribution theorem,— the chemical potential for free 
ions is given by 



In 



Pi 



In 



A^ 



(34) 



To obtain the corresponding expression for the chemical potential of ion pairs we use the 
Bethe approximation .^^iiSi^i In this case it yields 

(2p*/co) [1 - {2pl/c,)] 



Z2 



;i - Pi - 2p^)' 



(35) 



and, finally. 



Using the expression for the mass-action law and Eqs 



A« 



P.-H P^-^"'.-Pl^'' )+.n(g|. (36) 



, (jH^ and (jnni), we obtain 



Co 

ft = T 



1 - Co(pI)2exp 



27ra3 



1 



(37) 



This expression specifies the density of ion pairs in terms of densities of free charged particles. 

Because the ion pairs are neutral, in DHBj approximation they do not contribute into 
the electrostatic free energy. Thus, 



where the ideal contribution to free energy is given by^ 



(38) 



-f {pl^pl)=pl\npl-{l-pl- 2pl) ln(l -pl-2pl)- p\ In p\ - {p\ - 3) ln(l - 2p\lc^). (39) 
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The phase coexistence curve for sc lattice in DHBj approximation is shown in Fig. 4. The 
critical temperature is slightly lower, while the critical density is much larger than in the pure 
DH approximation. But it also has an unphysical banana-like shape, and this is because of 
the neglect of electrostatic interactions of ion pairs with free charged particles, as explained 
for continuum and lattice models of electrolytes.i2*i^ 



D. Dipole-Ion Interactions for Lattices with ( = 1 

Clearly unphysical phase behavior in DHBj approximation is due to the fact that the 
properties of ion pairs as dipoles are not taken into account. As was shown earlier for 
the realistic description of ionic fluids it is important to consider the ion pairs as dipoles that 
interact with residual free charged particles. These solvation effects eliminate the unphysical 
phase behavior and provide a better agreement between the calculated critical parameters 
and the parameters estimated from the computer simulations.i^ii^ 

Let us consider an arbitrary dipole particle that occupies 2 neighboring lattice sites. Let 
assume that the positive charge of the particle is at the site (0,0,0), while the negative 
charge is at (1, 0, 0). The linearized lattice Poisson-Boltzmann equation ((Tj) in this case can 
be written in the following form, 

A^(r) = K^ip{r) - Ao6{r) - Ai6{r - (1, 0, 0)). (40) 

Comparing with Eq. ()12j). we note that the last 2 terms in this expression reflect the fact 
that the particle is a dipole, consisting of two opposite charges. The Eq. (jin|) can be solved 
as before, 

a= /• yl„ + /lie'"!-''"'") 



k 



with coefficients Aq and Ai determined from the following conditions: 



A^{r = (0, 0, 0)) = -— ^ = /€V(r = (0, 0, 0)) - Aq, (42) 

UVq 



and 



Aifir = (1, 0, 0)) = = /tV(r = (1, 0, 0)) - A,. (43) 

These equations lead to the expression for the average electrostatic potential ipj^^ due to all 
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ions except the positive ion fixed at (0, 0, 0) and negative ion fixed at (1, 0, 0), 




(44) 



Dvo 6 [6-x2 + 4,P(^)J Dvo 6 



where we defined an auxihary function (^(a;^), 



(45) 



The contribution to the free energy density due to the interactions between dipoles and 
free ions can be calculated from 



where / {pi, pi) is given by Eq. (jSHl)- From these equations all thermodynamic properties 
can be obtained. 

The resulting phase coexistence curves and the critical parameters for sc, bcc and fee 
lattices estimated using the full DHBjDl approach are presented in Fig. 4 and Table 2, 
respectively. The predicted critical temperatures are decreasing as the number of nearest 
neighbors is going up, but they are still 25-30 % higher than the corresponding values for the 
continuum DHBjDl theoryi^ and for the Monte carlo simulations.— At the same time, 
the critical densities are approaching the predictions of continuum RPM calculations, while 
deviating significantly from the computer simulations estimates.— 

IV. SUMMARY AND CONCLUSIONS 

A theoretical investigation of ionic systems on the lattices with different discretizations 
is presented. A general lattice theory, based on the Debye-Hiickel approach, is developed for 
lattices with any discretization parameter. The theory is applied for explicit calculations of 
thermodynamic properties, phase behavior and critical parameters of several lattice models 
of electrolytes. 




(46) 







The overall free energy density of the system can be written as 



f = fipl,P*2) + f\pl)+f 



(47) 
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The simplest version of general lattice theory, the pure Debye-Hiickel theory that takes 
into account only the free ions, is utilized for obtaining the thermodynamic properties of sc, 
bcc and fee lattices for the discretization parameter C = I5 and for sc lattices with ( = -\/2, 
and 2. All considered cases exhibit the low-density gas-liquid phase transitions. As ex- 
pected, with the increase in the lattice discretization parameter ( the critical parameters are 
generally decreasing and approaching the values obtained in the continuum RPM treatment 
of ionic fiuids^^ that corresponds to C = 00 case. The behavior of the critical parameters 
is non-monotonic, especially for non- integer (, in agreement with observations from latest 
Monte Carlo computer simulations for large discretization parameters {( > 5).— 

For ionic fluids on the sc, bcc and fee lattices with ( = 1 a. more realistic theoretical 
description, that takes into account the creation of ion pairs and the ion-dipole interactions, 
is developed. The consideration of ion pairs as neutral species (the DHBj approach) slightly 
lowers the critical temperature and significantly increases the critical density, while the over- 
all phase coexistence curves have unphysical banana-like shapes. This unphysical behavior 
is cured when the ion pairs are treated as dipoles that interact with the residual free charged 
particles. The resulting critical temperatures decrease even more, while the values of critical 
densities are also lower. The critical temperatures for the lattice electrolytes are slightly 
higher than the continuum RPM values and the computer simulations estimates. The pre- 
dicted critical densities for the lattice models are essentially the same as in the corresponding 
continuum models of electrolytes. However, they are much lower than the critical densities 
from Monte Carlo simulations. This behavior is typical for lattice models. 

Our theoretical approach allows to calculate the thermodynamic properties of ionic sys- 
tems at low densities. However, at higher densities the sublattice ordering becomes very 
important in sc and bcc lattice models of electrolytes.'''^ The original lattice Debye-Hiickel 
approach^ was able to capture the ordering processes. It was predicted that the sublattice 
ordering would suppress the gas-liquid phase coexistence and the tricritical point would 
appear, in accord with the recent lattice computer simulations.-^ The predicted tricritical 
density was in agreement with the estimates from Monte Carlo computer simulations,"^^ 
while the critical temperature was overestimated. We plan to investigate the sublattice or- 
dering using the modified lattice DH method, developed in this paper, in a near future. In 
addition, our theoretical method might be used to investigate many other problems related 
to ionic fluids, such as the effect of additional short-range interactions, the size asymmetry 
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and charge asymmetry, and the multi-component mixtures. 
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Table I: Critical parameters for lattice electrolytes in the pure DH theory. HC corresponds to the 
predictions from the continuum RPM with hard-core interactions!^ 



Model 


C 




sc 


0.084377 


0.007978 


bcc 


0.083793 


0.005659 


fee 


0.082415 


0.004591 


HC 


0.061912 


0.004582 
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Table II: Critical parameters from the full DHBjDI theory for lattice and continuum electrolytes. 
For comparison, the estimates of Monte Carlo computer simulations are also shown. 



Model 


-'-c 


P*c 


sc 


0.0756 


0.0422 


bee 


0.0752 


0.0336 


fee 


0.0686 


0.0292 


continuun^i^ 


0.0522-0.0554 


0.0244-0.0259 


simulationo^i^^i^^ 


0.049 


0.06-0.08 
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Figure Captions: 



Fig. 1 A schematic picture of tlie lattice model of electrolytes with different discretizations. 
Large empty circles denote the positive particles, while large gray circles correspond to the 
negative ions. The ions are spheres with the diameter a(, where a is the distance between 
the nearest lattice sites. The position of the origin is marked by a cross. Small filled circles 
indicate the position of lattice sites in the exclusion zone of the central positive ion. 

Fig. 2 Phase coexistence curves predicted by pure DH theory for lattice electrolytes with 
C = 1: (a) sc, (b) bcc, (c) fee. The predictions from the continuum model of electrolytes 
with hard-core interactions^^ are shown in (d). 

Fig. 3 Critical parameters of lattice models of electrolytes for pure DH theory as a function 
of inverse lattice discretization parameters (: a) critical temperatures, b) critical densities. 

Fig. 4 Phase coexistence curves of ionic systems calculated for (a) sc lattice in the DHBj 
approximation; and using the full DHBjDI approach for (b) sc, (b) bcc, and (d) fee lattices. 
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